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We study, within the statistical hadronization model, the influence of narrow strangeness carrying 
baryon resonances (pentaquarks) on the understanding of particle production in relativistic heavy 
ion collisions. There is a great variation of expected yields as function of heavy ion collision energy 
due to rapidly evolving chemical conditions at particle chemical freeze-out. At relatively low collision 
energies, these new states lead to improvement of statistical hadronization fits. 
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Enhanced production of strange hadrons in relativistic 
heavy ion collisions is well established 0, • The avail- 
ability of a high abundance of strangeness favors pro- 
duction of strange hadron resonances, a topic of current 
intense experimental interest in the field of relativistic 
heavy ion collisions H EM IS 0, II The discovery by 
the NA49 collaboration (iq of a new S~ (1862) 1 = 3/2 
narrow T < 5 MeV resonance in their pp background, 
rather than in the AA foreground data at projectile en- 
ergy 158 GeV (^snn = 17.2 GeV) poses the question in 
which conditions one should look in heavy ion collisions 
for such new resonances. 

This newly discovered hadron resonance has, given 
the mass and charge, an exceedingly narrow width. 
This feature is common with 6 + (1540), another re- 
cently reported resonance 0, 0, 0, which de- 
cays into the channel with quark content uudds and 
1 = 0. This is believed to be the predicted [l5|. low- 
est mass, pentaquark state 0. The S*(1862) can be 
interpreted as its most massive isospin quartet member 
ssddu, ssudu, ssudd, ssuud with electrical charge vary- 
ing, respectively, from —2 to +1, in units of |e|. 

Appearance of these new resonances can have many 
consequences in the field of heavy ion collisions. We 
at first explore how the introduction into the fam- 
ily of hadronic particles of these two new resonances, 
9+(1540) and S*(1862), influence the results of statis- 
tical hadronization fit to relativistic heavy ion hadron 
production experimental results. We use the same data 
set as has been employed in Ref 0, El 0] and obtain 
predictions of how the relative abundances of these new 
resonant states vary as function of the heavy ion collision 
energy. 

Importantly, only the two already identified states with 
1 = 0, and / = 3/2 of the anti decuplet, which also 
includes the I = 1/2, and 1=1 states are of relevance in 
the study of the statistical hadronization fits. Thus, in 



our analysis, we do not depend on the unknown masses of 
/ = 1/2, and 1=1 states. However, the interpretation of 
the newly discovered narrow states as pentaquarks enters 
our considerations decisively. The pentaquark valance 
quark content enters the assigned chemical fugacities and 
phase space occupancies. The yield is proportional to the 
presumed spin degeneracy of the new states, taken to be 
two for a spin 1/2 anti decuplet. 

In our approach 0, 0, 0] , as in other recent work 
the chemical equilibrium and non-equilibrium is consid- 
ered. Accordingly, we allow quark pair phase space occu- 
pancies, for light quarks 7 g ^ 1, and/or strange quarks 
7 S / 1 [2(]]. Since we study at SPS the total particle 
multiplicities, and at RHIC the central yields which can 
be considered produced by rapidity-localized fireballs of 
matter, we require in our fits balance in the strange and 
antistrange quark content [2l| . 

There are two independent fit parameters when we as- 
sume complete chemical equilibrium, the chemical freeze- 
out temperature T and the light quark fugacity X q = 
V A u Ad = e Mi, /( 3T ). The baryochemical potential /ib 
is the physical parameter controlling baryon density. 
Strangeness conservation fixes the strange quark fugac- 
ity X s (equivalently, strangeness chemical potential, for 
more details see, e.g., @). Adding the possibility that the 
number of strange quark pairs is not in chemical equilib- 
rium, 7 S ^ 1, we have 3 parameters, and allowing also 
that light quark pair number is not in chemical equilib- 
rium, we have 4 parameters. These three alternatives 
will be coded as open triangles, open squares and filled 
squares, respectively, in all results we present graphically. 

We find that the new resonance + (154O) influences 
significantly the statistical hadronization fit to particle 
production at the lowest SPS energies. In a baryon rich 
environment the introduction into the fit of + (154O), 
a b = 1 baryon with 'wrong' strangeness influences the 
strangeness balance condition, and thus indirectly the 
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TABLE I: The chemical freeze-out statistical parameters found for nonequilibrium (left) and semi equilibrium (right) fits to 
SPS results. We show ^/sjvjv, the temperature T, light quark fugacity A 9 , strange quark fugacity A s , the quark occupancy 
parameters y q and 7s/7<j. Bottom line presents the statistical significance of the fit. The star (*) indicates for X s that it is 
a value resulting from strangeness conservation constraint. For 7, that there is an upper limit to which the value converged, 
7q < e m7r ^ T (on left), or that the value of 7 ? = 1 is set (on right). 



y/sNN [GeV J 


17.2 12.3 8.75 


17.2 12.3 8.75 


T [MeV] 


135 ±3 135 ±3 133 ±2 


157 ±4 156 ±4 154 ±3 


Xq 
Xs 


1.69(5) 1.98(6) 2.56(6) 
1.23* 1.27* 1.31* 


1.74(5) 2.03(7) 2.69(8) 
1.20* 1.24* 1.24* 


7g 


1.68* 1.68* 1.69* 
0.91(6) 0.83(4) 0.85(6) 


1* 1* 1* 
0.66(4) 0.60(4) 0.67(5) 


x7dof 


11.4/6 4.3/2 2.3/4 


23/7 8.9/3 4.0/5 



TABLE II: The chemical freeze-out statistical parameters 
found for nonequilibrium (left) and semi equilibrium (right) 
fits to RHIC results. We show ^/sjvjv, the temperature T, 
light quark fugacity X q , strange quark fugacity A s , the quark 
occupancy parameters 7, and Js/jq- Bottom line presents 
the statistical significance of the fit. The star (*) indicates for 
A s that it is a value resulting from strangeness conservation 
constraint. For 7, that there is an upper limit to which the 
value converged, 7, < e m7r ^ T (on left), or that the value of 
7q = 1 is set (on right). 



yJs NN [GeVj 


200 130 


200 130 


T [MeV] 


142 ±5 143 ±3 


159 ±6 159 ±2 


Xq 

X s 


1.051(9) 1.069(8) 
1.018* 1.023* 


1.052(9) 1.067(8) 
1.018* 1.023* 


7q 

Isllq 


1.62* 1.63* 
1.23(12) 1.32(5) 


1* 1* 
1.013(6) 1.13(4) 


X 7dof 


2.9/6 16.9/20 


4.6/7 32.7/21 
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FIG. 1: (Color online) x 2 /Aoi for statistical hadronization 
fits at SPS and RHIC: results are shown for 40, 80, 158A 
GeV Pb on stationary Pb target collisions and at RHIC for 
65+65 and 100+100A GeV Au-Au head on interactions. 



individual yields of all strange hadrons. This leads to a 
reduction in the statistical fit error for our hadronization 
study of the 40A GeV Pb-Pb reactions where we see a 
significant change in the relative yield of kaons and A. 
We also find changes in the details of the statistical fit 
parameters. In comparison to |17|. aside of the intro- 
duction of the new resonances, we also have harmonized 
our hadron decay table with those used by the Krakow 
group |22j . The improvement of the particle yield fit is 
both, a theoretical confirmation of the validity of the sta- 
tistical hadronization model of particle production, and 
its applicability at low SPS energies. 

We show how the fit error evolves in figure ^ which is 
also presented in the bottom lines of tables HI andlrfl along 
with the number of data points and resulting degrees of 
freedom. Considering the small number of degrees of 
freedom at SPS, we need x 2 /dof < 1 to have good signif- 
icance of the fit. The errors seen in figure Q are, for the 
chemical nonequilibrium case (filled squares), sufficiently 
small to allow us to conclude that the introduction of 



6 + (1540) assures that the statistical hadronization works 
well down to the lowest SPS energies. To compare with 
earlier results on x 2 /dof, obtained prior to the discovery 
of these new resonances, see Ref. [T3, figure 16. 

An interesting point, seen in figure^ is that the chem- 
ical equilibrium fit 7 S = 1, 7 9 = 1 is rendered unac- 
ceptable at all SPS energies in presence of the new reso- 
nances. The semi-equilibrium fit, which allows a varying 
strangeness saturation, but assumes light quark equilib- 
rium is generally resulting in twice as large x 2 compared 
to the full non-equilibrium approach. In a study of x 2 
profile as function of "f q we find a clear and strong mini- 
mum for 7 g — » <-y" lax = e -"W(2T)_ Acquisition by the fit 
of this limiting value implies that there is no fitting error 
in the "f q presented below. 

The chemical freeze-out parameters of the fits consid- 
ered play a very important role in predicting the (rela- 
tive) yield of hadronic particles, and this dependence is 
even stronger for many pentaquark states, due to their 
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FIG. 2: (Color online) Yield of 9+ (1540) in relativistic 
heavy ion collisions, based on statistical hadronization fit to 
hadronization parameters at SPS and RHIC 40, 80, 158A GeV 
Pb on stationary Pb target collisions and at RHIC for 65+65 
and 100+100A GeV Au-Au head on interactions. Relative 
yields with K a , A, and A(1520) are shown from bottom to 
top. 



unusual quantum numbers. These fit parameters for 
RHIC arc shown in table UH and for SPS in table [fl along 
with the freeze-out temperature. We note that for the full 
chemical non-equilibrium, the freeze-out temperature is 
found to be smaller than for semi-equilibrium case. This 
reduction is over-compensated in pentaquark yields by 
the significantly increased value of "/ q . 

We now consider the relative yields of the new reso- 
nances in figures [5] and [3] These yields vary strongly 
with collision energy for the case of + (154O) in figure[21 
but are rather constant in figure |31 Certainly our result 
differs greatly from expectations arising from an earlier 
study of the statistical model production of the + (154O) 
resonance [23| where the decisive variation of the particle 
yield with chemical potentials was not exp lored. More- 
over, the hadron yields, presented in |23|, did not in- 
clude the contributions from decay of short lived hadron 
resonances. We checked that the relative particle yields 
shown in piij for zero chemical potentials and varying 
temperature are mathematically correct, also as a cross 
check of our program. 

In figure [3 we show (from top to bottom) the relative 
yields 6+(1540)/A(1520), 0+ (1540) /A, 0+(154O)/^ s 
for chemical noncquilibrium (solid lines), semi- 
equilibrium ("f q = 1, dashed lines) and chemical 
equilibrium (dotted lines). The yields of A used here 



FIG. 3: (Color online) Relative yields H""(1862)/H" and 
E~(1786)/A are shown from bottom to top, see figure |U for 
further details. 



include 50% weak interaction cascade from S. 

The reason that the chemical nonequilibrium is lead- 
ing to greater than equilibrium yields is that the 
lower hadronization temperature is overcompensated 
by the chemical factors, e.g., 0+(154O)/A(152O) ~ 
l/2j q (X q /X s ) 2 ignoring the small mass difference. The 
factor 1/2 is due to the difference in spin degeneracy. The 
actually observed yield ratio 0+(154O)/A(152O) could be 
still greater since A(1520) is seen at 50% of the expected 
statistical hadronization yield in heavy ion collisions. In 
figure |21 we also recognize that the reason that there is 
such a significant impact at low SPS energies of + (154O) 
is that it is produced at the level of +10-20% of A in fits 
at 40A GeV. This is due to the large prevailing bary- 
ochemical density. Clearly, this is the environment in 
which one would want to study the properties of this new 
resonance in more detail. However, at all energies consid- 
ered, we find that + (154O) is more abundant compared 
to A(1520) and thus this new resonance could become an 
important probe of the hadronization dynamics. 

The observation of the pattern of relative yield of 
+ (154O), seen in figure |21 would firmly confirm the 4- 
quark, one anti quark content of this state. Namely, were 
for example the + (154O) another tri-quark baryon state, 
the yield ratio with (strange) baryons would be quite flat 
as function of collision energy. We further note that the 
absolute magnitude of the relative yield, seen in figure 
12 will be of help in establishing the degree of chemical 
equilibration. 

In figure |3 we show at the bottom the expected rel- 
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ative yield of the (1862) [ssddu] relative to H~ [ssd]. 
The S*(1862) adds at the percentile level to the yield of 
observed S and thus it is less influential in the statistical 
hadronization approach. The absence of variation of the 
relative yield with collision energy is due to cancellation 
of chemical factors. This relatively small relative yield at 
all collision energies here considered shows that indeed 
the pp environment, where it has been identified by the 
NA49 collaboration, is most suitable. The dotted lines, 
in figure |3 are visibly breaking the trend in some of the 
results, indicating that the large \ 2 chemical equilibrium 
fit generates unreliable statistical model parameters. 

We also show, in figure [3] on the top, the yield of the 
pentaquark state £(1776) [sddira] which for purpose of 
this study is assumed at the mass indicated. Again due 
to cancellation of key chemical factors in ratios shown in 
figure|3| both being proportional to -fq\d/X u , the ratio is 
flat (except for the failed fit chemical equilibrium results). 
Considering that ~ X u and 7 g ~ 1.6, the magnitude 



of relative yields seen in figure is primarily due to the 
hadron mass, and degeneracy. 

We have shown that inclusion of the pentaquark states 
in the study of particle production in heavy ion collisions 
improves the quality of our fits to experimental data. We 
find that 9 + (1540) state influences the low energy SPS 
particle yield fit results. It can be expected that it will be 
detectable, in particular at low heavy ion collision ener- 
gies, and thus should become a new probe of hadroniza- 
tion dynamics. The other pentaquark states will be hard 
to observe in heavy ion collisions. 
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